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Why to study quarkonia 

Quarkonia, i.e. bound states made of a heavy quark and 
a heavy antiquark (like charmonia, bottomonia, ...) are sys- 
tems where low-energy QCD may be studied in a system- 
atic way (e.g. one may address issues like large-order per- 
turbation theory, non-perturbative matrix elements, QCD 
vacuum, exotica, confinement, deconfinement, ... ). This 
is because M » p > E, where M is the heavy-quark 
mass, p the momentum transfer and E the binding energy 
of the bound state, and because M ^S> Aqcd, the scale of 
non-perturbative QCD. 

(1) M ^> p ^> E implies that quarkonia are non- 
relativistic and characterized by the hierarchy of 
scales typical of a non-relativistic bound state: 

M >p~ l/r~ Mv^>E~Mv 2 , (1) 

where r is the typical distance between the heavy 
quark and the heavy antiquark and v <C 1 is the typi- 
cal heavy-quark velocity. Systematic expansions in v 
may be implemented at the Lagrangian level by con- 
structing suitable effective field theories (EFTs): ex- 
panding QCD in p/M and E /M leads to NRQCD Q] ; 
expanding NRQCD in E r leads to pNRQCD Q. The 
hierarchy of non-relativistic scales makes the very dif- 
ference of quarkonia from heavy-light mesons, which 
are characterized by just two scales: M and Aqcd 0- 
For a review of non-relativistic effective field theories 
we refer to [4 |. Some historical background has been 
discussed in 0. 

(2) M 3> Aqcd implies a s (M) <C 1: phenomena hap- 
pening at the scale M may be treated perturbatively. 
We may further have small couplings if M v ^> Aqcd 
and Mv 2 ^> Aqcd, in which case a s (Mv) <C 1 and 
a s (Mv 2 ) <C 1 respectively. This is likely to hap- 
pen only for the lowest charmonium and bottomonium 
states (see Fig. [T}. 

It is precisely the rich structure of separated energy 
scales that makes quarkonium an ideal probe of confine- 
ment and deconfinement. The different quarkonium radii 
are differently sensitive to the Coulombic and confining in- 
teraction (see Fig. [2j». Hence, different quarkonia will dis- 
sociate in a medium at different temperatures, providing a 
thermometer for the plasma |7] [8) . 

All quarkonium scales get entangled in a typical ampli- 
tude involving a quarkonium observable, see Fig. [3] In 
particular, quarkonium annihilation and production happen 
at the scale M, quarkonium binding happens at the scale 
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Figure 1: a s running at one loop and a s (M j / ^v j / ^) and 

a s (M r{1S) v T{1S) ). 
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Figure 2: Static potential vs quarkonium radii taken from 

E). 

Mv, which is the typical momentum exchanged inside the 
bound state, while very low-energy gluons and light quarks 
(also called ultrasoft degrees of freedom) live long enough 
that a bound state has time to form and, therefore, are sen- 
sitive to the scale Mv 2 . Ultrasoft gluons are responsible 
for phenomena like the Lamb shift in QCD. 

Physics at the scale M 

Quarkonium annihilation and production happen at the 
scale M. The suitable EFT is NRQCD, which follows from 
QCD by integrating out the scale M, see Fig. |4] As a con- 
sequence, the effective Lagrangian is organized as an ex- 
pansion in 1/M and a s (M): 

£ N RQCD = 2^ jj^ X °n(^ Mv , Mv > -). 

n 

(2) 

where O n are the operators of NRQCD that live at the low- 
energy scales Mv and Mv 2 , /i is the NRQCD factorization 




Figure 3: Typical scales appearing in a quarkonium annihi- 
lation diagram. 



scale and c„ are the Wilson coefficients of the EFT that 
encode the contributions from the scale M and are non- 
analytic in M. 
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Figure 4: Matching of the diagram of Fig.[3]to NRQCD. 

The NRQCD factorization formula for the quarkonium 
annihilation width (into light hadrons or photons or lepton 
pairs) reads: 
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where do n is the dimension of the four-fermion operator 
£)4-fermion Recently, substantial progress has been made 
in the evaluation of the factorization formula at order v 7 
(9] [TO), in the lattice evaluation of the NRQCD matrix el- 
ements (iTIC^-fermion^ |TJJJ an( j j n me data of many 

hadronic and electromagnetic decays fl2l . As it was al- 
ready discussed in [12], the data are clearly sensitive to 
NLO corrections in the Wilson coefficients c n (and pre- 
sumably also to relativistic corrections). For an updated list 
of ratios of P-wave charmonium decay widths, see Tab.Q] 
The high precision of data and matrix elements has been 
recently exploited to provide a new determination of a s 
fromr T(ls) ^ 7Lh ./r T(15) ^ lh . IE): 



0.015 



a 8 (M r(ls) ) = 0.184+™ 



a s (M z ) = 0.119±™° l 
(4) 

Recent progress in quarkonium production in the frame- 
work of the NRQCD factorization have been discussed in 

ma. 

Physics at the scales Mv and Mv 2 

Quarkonium formation happens at the scale Mv. The 
suitable EFT is pNRQCD [2 1, which follows from NRQCD 
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Table 1: Comparison of decay width ratios of Xcj from 
Ifl3l (l.h. stands for light hadrons) with LO and NLO deter- 
minations (without corrections of relative order v 2 , m c = 
1.5 GeV and a s (2m c ) = 0.245). 
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Figure 5: Matching of the diagram of Fig. [4] to pNRQCD. 

by integrating out the scale Mv, see Fig. [5] As a conse- 
quence, the effective Lagrangian is organized as an expan- 
sion in 1/M and a s (M), inherited from NRQCD, and an 
expansion in r: 
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xV n , k {r, fi' , /i) r k x O k {n',Mv 2 ,...), (5) 

where Ok are the operators of pNRQCD that live at the 
low-energy scale Mv 2 , // is the pNRQCD factorization 
scale and V n ,k are the Wilson coefficients of the EFT that 
encode the contributions from the scale r and are non- 
analytic in r. Looking at the equations of motion of pN- 
RQCD, we may identify V n ,o with the 1/M n potentials 
that enter the Schrodinger equation and V Ul k^o with the 
couplings of the ultrasoft degrees of freedom, which pro- 
vide corrections to the Schrodinger equation. 

The static QCD spectrum without light quarks 

The spectrum of a static quark and a static antiquark has 
been studied in QCD without light quarks on the lattice, for 
instance in ifTBI . At short distances, it is well described by 
the Coulomb potential in the colour-singlet or in the colour- 
octet configurations: V s = — 4a s /3r and V = a s /6r. At 



large distances, the energies rise linearly with r. Higher ex- 
citations develop a mass gap of order Aqcd with respect to 
the lowest one. Reintroducing the heavy-quark mass M, 
the spectrum of the Mv 2 fluctuations around the lowest 
state constitutes the quarkonium spectrum while the spec- 
trum of the M v 2 fluctuations around the higher excitations 
constitutes the hybrid spectrum. 

Quarkonia 

The energy of the lowest excitation between a static 
quark and a static antiquark is the quarkonium static en- 
ergy. Quarkonia may be identified with the Mv 2 spectrum 
that differentiates once l/M corrections (first of all the ki- 
netic energy) are added to the effective Lagrangian, which 
reads 

£ p nrqcd = Jd 3 rS^ (ido-tL-V. + ^S 

+ultrasoftd.o.f., (6) 

where S is the colour-singlet quarkonium field. 

At short distances, the static potential is well described 
by perturbation theory up to NNNLL accuracy (for a recent 
analysis see lfT7]| ). Therefore the lowest-lying quarkonium 
states have a radius that is small enough for perturbation 
theory to apply (see Fig. [TJ). Higher-order corrections to 
the spectrum, masses and wave functions have been calcu- 
lated in IfTSl . Non-perturbative corrections are small and 
encoded in (local or non-local) condensates. Many param- 
eters and observables of the lowest quarkonium states have 
been calculated. They include: c and b masses at NNLO 
and partially at N 3 LO (for a review see [ 12 ], a more recent 
result is in [ 19 1); B c mass at NNLO |20|; B*, r\ c , r\\, masses 
at NLL quarkonium IP fine splittings at NLO El ; 
T(15), rjb electromagnetic decays at NNLL J23); T(15) 
and J/ip radiative decays at NLO l25lk T(15) — » 7776, 
J/ip — ► 7?7c at NNLO |24j: tt cross section at NNLO 
1 26 1; leading thermal effects on a quarkonium in medium: 
masses, widths, [27|, ... . For some recent reviews, we refer 
to H3|29). 

High-lying quarkonia are the Mv 2 fluctuations around 
the long-range tail of the potential. The long-range tail of 
the potential is not accessible by perturbation theory (see 
Fig. [TJ- However, the potential may still be expanded in 
l/M and each term of the expansion can be expressed in 
terms of field-strength insertions on a static Wilson loop 
1 30 1, which can be calculated on the lattice OTI . The re- 
sulting potential may be used in ©. The solution of the 
corresponding Schrodinger equation provides the quarko- 
nium masses and wave functions. A trivial example of ap- 
plication of this method is the mass of the h c . The lat- 
tice data show a vanishing long-range component of the 
spin-spin potential so that the potential appears to be en- 
tirely dominated by its short-range, delta like, part. This 
suggests that the 1 P\ state should be close to the centre 
of gravity of the 3 Pj system. Indeed, the measured mass 
of the h c by CLEO, M hc = 3524.4 ± 0.6 ± 0.4 MeV 



021 and E835, M hc = 3525.8 ± 0.2 ± 0.2 MeV l33l . 
is just on the top of the mass of the l 3 Pj centre of gravity: 
M c . . g .(l 3 P/) = 3525.36 ±0.2 ±0.2 MeV. 

Gluonic excitations of quarkonia 

Many states, built on each of the hybrid potentials, are 
expected. These states typically develop a width also with- 
out including light quarks, since they may decay into lower 
states, e.g. like hybrid — > glueball + quark-antiquark. 

One possible candidate for such a state is the y(4260). 
The F(4260) has been discovered by BABAR in the ra- 
diative return process e + e~ — * 77r + 7r~ J/ip with mass 
M Y = 4259 ± 8±| MeV and width Y = 88 ± 23^ MeV 
(34J, and seen in the same process by BELLE with mass 
M Y = 4247 ± 12±\l MeV and width Y — 108 ± 19 ± 
10 MeV E3 and by CLEO with mass M Y = 4284+^ ± 
4 MeV and width Y = 73±|§ ± 5 MeV (36). CLEO 
has also confirmed the existence of an enhancement in 
the h + tt~ J/ip cross section at 4260 MeV in a measure- 
ment of direct e + e~ annihilation at s/s = 4040, 4160 
and 4260 MeV 071- The V(4260) J PC quantum num- 
bers are 1~. BABAR measures B(Y -> DD)/B(Y -> 
J/-07T+7T-) < 1.0 (« 500 for V(3770), which suggests 
an exotic interpretation for the V(4260)) [38], moreover 
B{Y_ -> D*D)/B{Y -> J/i/m+ic-) < 34 and B(Y -> 
D*D*)/B(Y -> J/tPit+tt-) < 40 BJJ. 
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Figure 6: The hybrid static potential II„ at short and in- 
termediate distances, ro ~ 0.5 fm. The points are the lat- 
tice data from [16], the continuous line is the fitting curve 
En u ro = constant + 0.11 r Q /r + 0.24 (r/r ) 2 . 

Many interpretations have been proposed for the 
F(4260), one of this is that the V(4260) is a charmo- 
nium hybrid igOJ STJ 02) . If the V(4260) is interpreted 
as a charmonium hybrid, one may rely on the heavy-quark 
expansion and on lattice calculations to study its proper- 
ties. Decays into D^D^ should be suppressed, since 
they are forbidden at leading order in the heavy-quark ex- 
pansion fiD . This is in agreement with the upper limit on 
Y — > DD reported by BABAR. The quantum numbers of 
the V(4260) are consistent with those of a pseudoscalar 
^ fluctuation \<fi) belonging to the family of Mv 2 fluctu- 
ations around the gluonic excitation between a static quark 
and a static antiquark with quantum numbers l" 1 , also 



known as II U , 

\y) = |n u ) ® (7) 

It is suggestive that, according to lattice calculations [16|, 
n„ is the lowest gluonic excitation between a static quark 
and a static antiquark above the quark-antiquark colour sin- 
glet. \(p) may be obtained as the solution of the Schrodinger 
equation whose potential is the static energy of n„. Fitting 
the static energy of U u at short and intermediate distances, 
one gets -En„?"o = constant + 0.11 ro/r + 0.24 (r/ro) 2 , 
see Fig. [g] Solving the corresponding Schrodinger equa- 
tion, one gets My = 2 x 1.48 + 0.87 + 0.53 = 4.36 GeV, 
where 1.48 GeV is the charm mass in the RS scheme [43] 
and 0.87 GeV is the gluelump mass in the same scheme 

The QCD spectrum with light quarks 

Adding light quarks changes the heavy quark-antiquark 
spectrum in the following way (see J. Soto in l45lD . 

(1) We still have states made just of heavy quarks and glu- 
ons. They may develop a width because of the decay 
through pion emission. If new states made with heavy 
and light quarks develop a mass gap of order Aqcd 
with respect to the former ones, then these new states 
may be absorbed into the definition of the potentials 
or of the (local or non-local) condensates [46 1. 

(2) In addition, new states built using the light quark 
quantum numbers may form. Possible states made of 
two heavy and light quarks include states built on the 
pair of heavy-light mesons (DD, BB, ...), molecular 
states [47 1, molecular states made of the usual quarko- 
nium states, built on the static potential, and light 
hadrons (hadro-quarkonium) [48], pairs of heavy-light 
baryons ||49ll , tetraquark states |[50l and likely many 
other states. 

Clear evidence for four-quark states may be provided by 
a charged resonance, like the Z+(4430), Z+(4050) and 
Z^~(4250) signals, detected by BELLE (but not confirmed 
by BABAR), possibly are. See EJUS and the panel dis- 
cussion at this conference. 

There is accumulating evidence, although not yet con- 
clusive, that the A (3872) may be a four quark state. 
The state A (3872) has been discovered by BELLE in 
B ± -> K ± X -> K ± tt + tt~ J/ip with M x = 3872.0 + 
0.6 ± 0.5 MeV (53), and confirmed by BABAR El that 
measures M x = 3871.4 ± 0.6 ± 0.1 MeV in B+ -> 
K+ir+n-J/xJj and M x = 3868.7 ± 1.5 ± 0.4 MeV in 
B° — > K°n + iT~ J/ijj Il55l . The state has also been seen 
at the Tevatron in pp — + A — > n + ir~J/ip by CDF with a 
mass M x = 3871.3 ± 0.7 + 0.4 MeV Eg) and by DO with 
a mass M x = 3871.8 ± 3.1 ± 3.0 MeV E3- BELLE 
has an upper limit on the width: V < 2.3 MeV, while 
BABAR finds T = 3.0±^ ± 0.9 MeV El. The A(3872) 
has been detected in different decay modes, the decay 
X — > D°D°tt being likely the dominant one: B(X — > 



D°D 7r°)/B(X -» tt+tt-J/V^) = 9.4±|;| E3. One 
should notice that BELLE finds a threshold enhancement 
peak in the D°D°ir° invariant mass at 3875.4 ± 0.7±l'% 
MeV, which is about 2 a larger than the world-average 
mass of the A(3872). The decay mode X -> 7 J/ip IgO) 
implies that the A (3872) has positive charge conjugation. 
Analyses of angular distributions performed by BELLE 
[61 1 and CDF |62| favor a spin parity assignment 1 + . The 
ratio B{X -> tt+tt"^ J/ip)/B(X -► ir+ir-J/iP) = 
1.0 + 0.4 + 0.3 measured by BELLE (60) suggests that 
the A(3872) is a mixture of isospin 1=1 and 1 = 
states. The substantial 1 = 1 component requires that the 
A (3872) contains uuldd pairs in addition to hidden charm, 
which thus qualifies it as a four-quark state ll63l . Hence, 
most recently, the majority of theoretical studies has ana- 
lyzed the A (3872) as a four-quark state with J PC quantum 
numbers 1 ++ . See Fig. [7] 




Figure 7: ccqq 4-quark state. 

Three quark-pair configurations are possible for a four- 
quark state of the type ccqq (q stands for a generic light 
quark). All of them have been exploited in the litera- 
ture. However, the resulting models are not equivalent, be- 
cause different dynamics are attributed to different confi- 
gurations. 

Clearly, for states made of heavy quark-antiquark pairs 
and light quarks, it would be useful to have the spectrum of 
tetraquark potentials, like the one that we have for the glu- 
onic excitations and that we discussed above. It would al- 
low us to build a plethora of states on each of the potentials, 
many of them developing a width due to decays through 
pion (or other light hadron) emission. Diquarks have been 
recently investigated on the lattice ||64) . In the lack of fur- 
ther theoretical input from QCD, many tetraquark studies 
rely on phenomenological models of the tetraquark inter- 
action unless some special hierarchy of dynamical scales 
may be further exploited on the top of the non-relativistic 
and perturbative expansions discussed so far. 

In l65l . see Fig. [8] it is assumed that A ~ 
(cc) | = i(<7<7)! = i, i.e. that the dominant Fock-space com- 
ponent contains a cc pair and a qq pair in a colour-octet 
configuration with spin 1 . Calculations have been based on 
a phenomenological interaction Hamiltonian. 

In [66 1, see Fig. [9] it is assumed that A ~ 
(cg)| =1 (cg)! =0 + (cg)| =0 (cg)! =1 . Here, the clustering 
of quark pairs in tightly bound colour-triplet diquarks is 
not induced by a scale separation as it would happen in 
baryons made of two heavy quarks l67ll . but is a dynami- 
cal assumption of the model. In particular, the model pre- 




Figure 8: ccqq 4-quark state in the configuration of [65 1. 



diets the existence of two neutral states made of cucu (X u ) 
and eded (A^) and of two charged ones. The two reso- 
nances discovered by BELLE and BABAR, the first decay- 
ing in J/ip7r + iT~ and the second preferably in D°D 7r° 
have been suggested as possible candidates for the Xd and 




Figure 9: ccqq 4-quark state in the configuration of |66 1. 

In (69] ED), see Fig. [TDJ it is assumed that X ~ 

(c?)ko(?g)ki + (c?)ki(95)ko ~ D ° D*° + D*°D°, 
i.e. that the dominant Fock-space component of the 
A (3872) is aD°D*° and D* D° molecule; small short- 
range components of the type (cc)g =1 (qq)g =1 ~ J/ip p, to 
are included as well. Predictions depend on the adopted 
phenomenological Hamiltonian, which typically contains, 
in the short range (~ 1/ AqcdX potential-type interactions 
among the quarks and, in the long range (~ 1/m,), the 
one-pion exchange. The prediction T(X — > ir + ir~ J ftp) « 
T(X — > ir + ir~ir° J/tp) made in ll70l turned out to be con- 
sistent with the BELLE result |60|. However, another pre- 
diction, T(A -> tt+tt-J/V') w 20T(A -> D°D°ir ), 
is two orders of magnitude far from the data ||59ll . Not 
necessarily this points to a failure of the molecular model, 
but possibly to a smaller J /tp p component in the A (3872) 
Fock space. 




Figure 10: ccqq 4-quark state in the configuration of |69 

ED. 



In ||7T1 l72l l73l . see Fig. [TT] it is assumed not only that 
the A(3872) is a D° D*° and D° D* molecule, but also 
that it is loosely bound, i.e. that the following hierarchy of 
scales is realized: Aqcd S> m, ^> rn^/Mo w 10 MeV 
> -^binding- Indeed, the binding energy, ^binding, which 
may be estimated from Mj — (M D , o +M D o ), is very close 
to zero, i.e. much smaller than the natural scale m^./MD . 
This is also the case when using a recent CLEO determi- 
nation of the D mass, M Da = 1864.847 ± 0.150 ± 0.095 
MeV [ 74 1 . The main uncertainty comes from the A (3872) 
mass. Systems with a short-range interaction and a large 
scattering length have universal properties that may be ex- 
ploited: in particular, production and decay amplitudes fac- 
torize in a short-range and a long-range part, where the lat- 
ter depends only on one single parameter, the scattering 
length. The long-range molecular D° D* and D° D* 
components of the A (3872) should be responsible for the 
A (3872) decaying into D°D°tt°. For a recent analysis of 
the BELLE data about the D°D°tt° final state enhance- 
ment and the molecular picture we refer to f75l . For discus- 
sion about the evaluation of the A (3872) production cross 
section at the Tevatron inside the molecular model we refer 




Figure 1 1 : ccqq 4-quark state in the configuration of ||7T1 

Eats). 

Coupled channels 

An important (and yet unsolved) problem is how all the 
different kind of states (with and without light quarks) in- 
teract with each other. A systematic treatment does not ex- 
ist so far. For the coupling with two-meson states, most 
of the existing analyses rely on two models, which are now 
more than 30 years old: the Cornell coupled-channel model 
ll78l and the 3 Po model [79]. Steps towards a lattice based 
approach have been undertaken recently ll80l IHTI and may, 
in perspective, provide an alternative, QCD based, treat- 
ment. 

Conclusions 

Our understanding of how a (effective field) theory of 
quarkonium should look like has dramatically increased 
over the last decade. 

For states below threshold such a theory exists and al- 
lows a systematic study of the quarkonium lowest reso- 
nances. Even precision physics is possible. Higher reso- 
nances may need to be supplemented by lattice data. High- 
quality lattice data have become available in the last years 



for some crucial quantities (e.g. potentials, decay matrix 
elements, ...). 

For states above threshold, the picture appears much 
more uncertain. Many degrees of freedom seem to be 
present, and the absence of a clear systematics is an ob- 
stacle to an universal picture. Most likely, systematic de- 
scriptions will be found that suite only specific families of 
states, the near-threshold molecular states providing an ex- 
ample. 
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